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Abstract. We report on two recent calculations of QCD corrections to neutralino annihilation cross 
sections: (1) the next-to-leading order corrections to %% — > gg, and (2) the contribution to the cross 
■ section for XX ~ * QQg arising from interference between the tree-level and loop-induced processes. 

> 

o i 

The presence of non-baryonic dark matter in the universe is compelling evidence for 
physics beyond the Standard Model. Supersymmetry (SUSY) provides an especially 
attractive explanation with the lightest supersymmetric particle, usually a neutralino X, 
as the dark matter candidate. In this talk we review two recent calculations yj,|2j] of QCD 
\ corrections to neutralino annihilation via internal squark exchange. 

The tree-level neutralino annihilation process XX ~~ * qq y i a exchange of a squark of 
mass Mq can be reduced to an effective vertex described by a dimension-six operator 



O 



X 



suppressed by M 2 



J? = (c/M%)0 6 , ^=(OT)(?fM). (1) 



If the annihilation occurs at rest, the Majorana neutralino spinors reduce to a spin-singlet 
(D \ combination, XY^iYsX ~ > KuXijs/ V2m x )x, where K is the center-of-mass momentum 

of the two-neutralino system. The operator G§ can then be written as the divergence of 
the axial-vector current, 



X (iy 5 /V2m x )x] [d^qfysq)} • (2) 

In the massless quark limit, m q = 0, the axial vector current is conserved at tree level, 
du(qY^Y5 ( l) = 0, and all tree-level dimension-six amplitudes vanish in this limit. This is 
the well known partially-conserved axial current (PCAC) condition. 
The quark mass suppression can be lifted in two ways: 

1 . By going beyond leading order in a s to include the correction to the dimension-six 
operator that involves the anomalous triangle diagram. This is XX ~ > 88- This pro- 
cess first appears at order a 2 , and thus suffers from a large QCD scale uncertainty. 
We compute the next-to-leading order (NLO) QCD corrections to this process 10]. 

2. By going to dimension-eight by including hard gluon radiation. This is XX ~^ IVS- 
The tree-level diagram interferes with XX ~^ 88* wim m e off-shell gluon splitting 
into a qq pair, yielding a contribution to the cross section of order aj, the same 
order as the leading XX ~^ 88 cross section. We compute this interference term J2]. 



XX ~^ a t NLO: The leading-order (LO) amplitude for XX ~~ * 88 m me massless 
quark limit is given by 
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where the sum runs over the five light (massless) quark flavors c( and we neglect terms 
subleading in m|/M~. The %4^rl couplings g r ^ for right and left quark helicities, 
respectively, are 



lr = -y/lNug'Q, 81 = -VlNn8'(T 3 -Q) + V2N ng T 3 . 



(4) 



Here T 3 is the quark isospin, Q is the quark electric charge, g and g' are the weak 
couplings under SU(2)l and U(l)y, respectively, and iVii and N\2 are the bino and 
wino components of the neutralino as defined in Ref. |0]. The matrix element in Eq. © 
results in an annihilation cross section for XX ~* 88 °f 
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For 6 GeV < m x < 110 GeV, this approximate formula is within 10% of the result 
including the full quark mass dependence LI]- 

In order to compute the NLO QCD corrections to XX ~^ 88> we exploit the anomaly 
equation to relate this process to pseudoscalar decay to gluons. The conservation of the 
axial vector current is violated by the quark mass and by the anomalous triangle diagram, 
in the form J4j] 

M^Ysq) =2m q qiy 5 q + ^-Gf v G^\ (6) 



with \G^ V = £f 1V apG a P denoting the dual color field strength tensor. 

The amplitude for XX ~^ 88 ls obtained by setting m q = and neglecting the first 
term on the right hand side of Eq. ©. The gluonic decay amplitude of a fundamental 
pseudoscalar, A — > gg, is also related to the anomaly equation. This decay proceeds 
through a quark loop; if the Yukawa coupling of this quark to A is proportional to the 
quark mass, the decay is dominated by the contribution of the heaviest quark in the loop. 
In the heavy quark limit, mg 3> m^, we can neglect the divergence term on the left-hand 
side of Eq. ©. Thus we see that two seemingly different processes, A — > gg through a 
heavy quark loop and XX ~^ 88 through a light quark loop, are related by the anomaly 
equation to the same gluonic operator. The Adler-Bardeen theorem ||5|] guarantees that 
Eq. © is valid to all orders in a s . We take advantage of this to obtain the NLO QCD 
corrections to XX ~^ 88 from the known results for A — * gg at NLO. 

In the heavy top quark limit, for which our anomaly relation is valid, the NLO QCD 
corrections to the A — ► gg partial width are given by a multiplicative factor times the LO 
decay rate [6]. Invoking the Adler-Bardeen theorem, the NLO cross section for XX ~^ 88 



in the zero-velocity limit becomes Jl|] 



V re lONLo(ZZ ~* 8S) = VrelOdim6 
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where /i is the renormalization scale. The strong coupling is evaluated based on Nf- 
flavor running at the scale \l where it appears both explicitly in the correction term and 
within the leading-order cross section v re iCdim6- The integer Nf counts the number of 
quark flavors in the gluon splitting, with Nf = 5 for m^, <C m x <C m t . Taking \x = 2m x , 
N f = 5, and m x = 100 GeV, we obtain v re ia N Lo(ZZ -> gg) — v re iOdim6[l +0.62]. The 
NLO corrections reduce the residual renormalization scale dependence of the %% ~^ 88 
annihilation cross section from ±16% to ±9% when we vary jU by a factor of two in 
either direction Jl|]. 

XX —> QQg interference term: To leading order in m| /M~, it is straightforward to ex- 
tend the loop amplitude for XX ~^ 88 t0 include one off- shell gluon splitting into a 
quark-antiquark pair: 
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where q\, q 2 , and q?, are the momenta of the final-state quark, antiquark, and gluon, 
respectively, and the sum runs over the five light (massless) quarks in the loop. 
The matrix element for the tree-level process XX ~^ 1^8 1S 
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where we neglect terms subleading in m 2 x /M~. Only the vectorlike part of ^iree will 
interfere with the dimension-six amplitude. Summing over final-state polarizations and 
colors and averaging over initial-state polarizations, we obtain the interference term yO, 
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where we have summed over the five light internal quarks q 1 in ^# sp ii t of Eq. © and over 
the five light external quarks q. Integrating over the phase space, we find the contribution 
to the cross section from the interference term J2], 
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The form of this interference term is rather similar to that of the dimension-six cross 
section given in Eq. ©. In the special case of degenerate right- and left-handed squarks, 
Mq R = Mq h = Mg, we find that the interference term is related to the dimension-six 
cross section by a multiplicative constant: v re iOi n t = [— 2m^/3M|]v r eiC7dim6- Combining 
this result with the NLO corrections to XX ~ > gg w i tn Nf = 5 light flavors and a 
renormalization scale ji = 2m x , we obtain 



v K \<J = v re ia dim6 
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For m% ~ 100 GeV the NLO correction is roughly 60%; for, e.g., Mq ~ 2m x , the 
interference term cancels off about one quarter of the NLO correction. We note that 
the expansion of the squark propagator in XX — ► gg yields correction terms starting at 
order /M~ in the square brackets in Eq. (fT2"l) : these are beyond our present interest. 

Outlook: In the region of SUSY parameter space where virtual squark exchange gives 
a sizable contribution to neutralino annihilation, the neutralinos are "p-wave annihila- 
tors" during freeze-out in the early universe; i.e., they annihilate through the velocity- 
dependent XX ~^ <l<i process, which is not suppressed by the quark mass. In this situ- 
ation, the process XX ~> gg is subdominant and our QCD corrections constitute only 
about 1 — 2% of the dominant cross section. They are thus of little importance for com- 
puting the relic neutralino abundance. However, for neutralino annihilation at the present 
time the relative neutralino velocity is much lower, leading to a much smaller tree-level 
XX ~ > cross section. In this situation, the XX ~> gg cross section can be as large 
or larger than XX ~> <!<!■> so tnat tne QCD corrections computed here could have a sig- 
nificant impact on the computation of gamma ray and neutrino fluxes from neutralino 
annihilation in the galatic halo and inside the sun, respectively. In [2] we provided the 
spin-averaged squares of matrix elements for the dimension-six, dimension-eight, and 
interference term contributions to XX ~^ QQg m terms of the final-state particle momenta 
to facilitate their implementation into Monte-Carlo generators for the computation of 
indirect dark matter detection rates. 
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